arXiv: 1503.0891 lv2 [gr-qc] 29 Aug 2015 


Variations on an aethereal theme 


Ted Jacobson* and Antony J. Speranza^ 

Maryland Center for Fundamental Physics, University of Maryland, College Park, MD 20742 and 
Perimeter Institute for Theoretical Physics, 31 Caroline Street North, Waterloo, ON N2L 2Y5, Canada 

(Dated: June 16, 2015) 

We consider a class of Lorentz-violating theories of gravity involving a timelike unit vector field (the aether) 
coupled to a metric, two examples being Einstein-aether theory and Hofava gravity. The action always includes 
the Ricci scalar of the metric and the invariants quadratic in covariant derivatives of the aether, but the theories 
differ in how the aether is constructed from other fields, and whether those fields are varied in the action. Fields 
that are not varied define background structures breaking diffeomorphsim invariance, including threadings, fo- 
lations, and clocks, which generally produce novel degrees of freedom arising from the violation of what would 
otherwise be initial value constraints. The principal aims of this paper are to survey the nature of the theories 
that arise, and to understand the consequences of breaking diffeomorphism invariance in this setting. In a com¬ 
panion paper [arXiv: 1504.03305], we address some of the phenomenology of the “ponderable aether” case in 
which the presence of a background clock endows the aether with a variable internal energy density that behaves 
in some respects like dark matter. 


I. INTRODUCTION AND SUMMARY 

Longstanding puzzles of cosmology and quantum gravity 
have led some to question the fundamental assumption of gen¬ 
eral relativity, that the spacetime manifold has no structure 
other than that determined by the metric. In particular, the 
cosmological constant problem, dark energy, dark matter, the 
trans-Planckian puzzle, the need for a UV completion of gen¬ 
eral relativity, the problem of time and the interpretation of 
quantum cosmology have motivated exploration of modified 
gravity theories with vacuum structure violating local Lorentz 
boost symmetry. If exact rotational symmetry is preserved, a 
Lorentz violating vacuum structure selects a preferred time¬ 
like direction at each spacetime point. The integral curves of 
this field of directions may be thought of as the flow of an 
“aether fluid.” The 4-velocity of the aether is the unit time¬ 
like vector held tangent to this flow. 

In constructing a theory with such an aether, one must de¬ 
cide whether the aether is to be treated as dynamical, i.e. var¬ 
ied in the action principle, or instead as background. If the 
aether is dynamical, then one must further specify how it is 
constructed in terms of the fields that are varied in the action. 
Actually it turns out that the distinction between varied and 
not varied fields is not so clear cut: the equations of motion 
for scalar fields are often a consequence of the equations of 
motion of the other fields. Such scalars, and the structures 
they define, can therefore be regarded as “background” struc¬ 
ture, even though they might also be varied in the action. What 
is important for the physics, however, is not how we refer to 
them, but how these choices affect the degrees of freedom and 
behavior of the theories. The purpose of this paper is to ex¬ 
amine this question for a variety of related aether theories. 

It is natural to assume that the aether-metric dynamics is 
governed by an (effective) action involving the metric, its cur¬ 
vature, the aether, and covariant derivatives of the aether and 
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the curvature. Before beginning with the detailed analysis, 
we would like to point out that the theory would generally 
be dynamically overconstrained, i.e. “inconsistent,” if deriva¬ 
tives of the aether were not included in the action. Suppose for 
example the action S\gab: </>, u°‘] is a scalar constructed from 
the metric gab, a scalar field </>, and a unit vector field u“. 
Suppose further that the aether enters the Lagrangian density 
only via the coupling The variations of 

must be orthogonal to to preserve the unit condition, and 
these impose the equation of motion {SS/Su°'){S^ — u°'Ub) = 
y/^{u^(j)^rn)[(t>,b “ {u'^4‘,n)ub] = 0. This extremely restric¬ 
tive condition requires that either (p is constant along the flow 
lines of or the flow of is hypersurface orthogonal and (p 
is constant on the orthogonal hypersurfaces. This eliminates 
virtually all of the solutions to the scalar equation of motion. 
Moreover, even if we choose to not impose the aether equa¬ 
tion of motion, the scalar held is still overconstrained since, as 
shown below, the other equations of motion imply that the Lie 
derivative £u[{5S/5u°‘){5^ — u°‘Ub)] vanishes. Although a 
weaker condition, this still eliminates almost all scalar held 
solutions. The situation is quite different, however, if the 
action includes terms quadratic in aether derivatives. Then 
5S/5u°‘ includes second derivatives of the aether, so instead 
of overconstraining the scalar held the extra conditions can be 
propagation equations for the aether. 

At lowest order in a derivative expansion, the most general 
action for the metric and aether is given (up to the integral of 
a total divergence) by ' 

S\gab,U {c,9ab)] — 

^ / d'^Xy/^{R+ ^0“^ + CaCt'^ + CojUJ^ + CaO^), 

IbTTGo J 3 

( 1 ) 


^ We use the metric signature (H-). Abstract indices are denoted by 

Latin letters, spacetime coordinate indices by Greek letters, and comma 
and semicolon before an index denote partial and covariant derivative re¬ 
spectively. Quantities with density weight 1 are written in caligraphic font, 
£, C, or carry a tilde, unless they involve the metric determinant or are 
written explicitly as a variational derivative. 
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where C denotes collectively independently varied fields used 
in the construction of the aether gab)- The terms in the 
integrand are the Ricci scalar R of the metric, and the expan¬ 
sion 9 = Wau‘^, shear aab = V ^^Ub) + ^ {VcU^)hab - ui^aO-b), 
twist ujab = — itfaUb] and acceleration ab = u°‘y aUb of 

the aether flow (here hab = UaUb — gab is the spatial metric). 
In this paper we examine how variations on the construction 
of rt“(C, gab) affect the resulting theory. 

An important distinction is whether or not the aether vec¬ 
tor is necessarily hypersurface-orthogonal, and hence non¬ 
twisting. The minimal structure required to determine a twist- 
free aether is a foliation by spacelike hypersurfaces, whereas 
the minimal structure required to determine a twisting aether 
is a timelike congruence of curves, i.e. a threading [1, 2]. A 
well-known example of a foliation theory results if is con¬ 
structed from a scalar field T as 

u'^{T,gab) = g'^Xb/\dn ( 2 ) 

where |dr| = {g°‘'^TaTb)^^'^ is the norm of the gradient. By 
construction, the aether (2) is orthogonal to the constant T 
surfaces. This is khronometric theory [3-5], a.k.a. the infrared 
limit of the nonprojectable version of Hofava-Lifshitz gravity 
[6]. Among threading theories, a well-known example results 
if u°‘ is constructed from an independent vector field as 

= (3) 

where |A| = is the norm of A“. This is 

Einstein-aether theory [7, 8], written in the form given in ref. 

[9] . In (3) the unit constraint on u°‘ holds by construction, 
rather than being imposed via a Lagrange multiplier term as is 
more commonly done. 

In the constructions we study, the foliation or threading is 
in most cases described covariantly using scalar fields which 
enter the construction of through their gradients. Since 
the action (1) involves first derivatives of m“, it involves sec¬ 
ond derivatives of the scalars, raising the concern that the re¬ 
sulting theory might suffer from an Ostrogradski instability 

[10] . However, diffeomorphism invariance implies that the 
Lagrangian is degenerate, so that there is a possibility that the 
instability is absent. In fact, as we explain in section II, for all 
theories we consider, the scalar field equations are redundant 
with the other field equations. Therefore the scalars need not 
be varied in the action principle, so we may fix their values by 
a choice of coordinates at the level of the action. The gauge- 
fixed action is no longer invariant under the full group of dif- 
feomorphisms, yet it gives a theory that is equivalent to the 
one defined by the original diffeomorphism invariant action.^ 
There exists a coordinate gauge choice for which the scalars’ 
gradients have constant components, and in such a gauge the 
action is only first order in derivatives of the remaining dy¬ 
namical variables. (For example, in the khronometric theory 


^ It is notoriously challenging to define what it means for a theory itself —as 
opposed to its formulation—to be diffeomorphsim invariant [11, 12]. 


with one dynamical scalar T, in (2) is first order in deriva¬ 
tives, but with the gauge choice = T its components be¬ 
come y /which contain no derivatives.) We 

conclude from this that there is no Ostrogradski instability in 
any gauge. 

In Einstein-aether theory, the threading is determined by a 
“line field,” i.e. a vector field modulo local scaling. We call 
this the dynamical aether theory, since it arises from a dy¬ 
namical field A°‘ that appears with only first derivatives in the 
action. In subsection III A, we consider a related theory where 
the threading is determined by three scalar fields, that can be 
fixed as background structures as explained above. We show 
that t\\\s fixed threading theory is equivalent to Einstein-aether 
theory except that it admits violation of the spatial initial value 
constraints. The constraint violation is characterized by a spa¬ 
tial covector density that is preserved along the aether flow, 
and does not affect the energy-momentum tensor. 

Subsection IIIB considers a different theory, in which the 
aether threading is determined by a line field as in Einstein- 
aether theory, and there is an additional scalar field which 
determines a preferred clock constrained to measure proper 
time along the threads. Thm fixed clock theory is equivalent 
to Einstein-aether theory except that it admits violation of a 
single initial value constraint per spatial point. The constraint 
violation is characterized by a scalar density that is preserved 
along the aether flow, and appears in the energy-momentum 
tensor like a rest mass density of the aether fluid. We call 
this a ponderable aether, invoking the 19th century adjective 
that was used to distinguish ordinary matter from aether. Fi¬ 
nally in subsection III C we consider the fixed aether theory, 
containing both a fixed threading and a fixed clock. The four 
scalars that describe the background structure completely de¬ 
termine the aether vector, and hence this theory is equivalent 
to Einstein-aether theory with the vector taken to be non- 
dynamical. 

In section IV, we consider foliation-type theories. After re¬ 
viewing the khronometric theory, we consider in section IV B 
introducing an independent clock field. As before this leads 
to a violated constraint as well as an additional component in 
the stress tensor. Finally, we show that when the clock field 
is constrained to coincide with the preferred foliation, the re¬ 
sulting theory is projectable Hofava gravity, and we again find 
a violated constraint. This constraint violation was already 
studied in ref. [13], which referred to it as “dark matter as an 
integration constant.” Unlike before, however, the violation in 
this case is not preserved under flows of for a generic form 
of the action (1). 


II. VARIATIONAL PRELIMINARIES 

We begin by establishing the notation and several key re¬ 
sults that are used throughout this paper. We define the follow¬ 
ing tensor densities, resulting from variations of the Einstein- 
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aether action (1), 


£ab 


£c 


6S 

§gab 


Ta 


SS 

§gab 


Su'^ 



( 4 ) 

( 5 ) 


The object to the right of the vertical line in these expressions 
indicates which tensors are held fixed when computing the 
variation, i.e. £ab and J^ab differ in that contravariant is 
held fixed in the former, while covariant Uc is fixed in the 
latter. They are related by 


^ab — £ab T ^(a^b)- (^) 

As discussed in the introduction, we are interested in cases 
where is constructed from gab and other fields collectively 
denoted (. The Einstein equation results from varying the ac¬ 
tion holding ( fixed, and hence receives a contribution from 
the explicit metrics appearing in the action, as well as from 
the metric dependence of In all cases considered, de¬ 
pends algebraically on the metric, and the Einstein equation 
takes the form 

£ab+£e^=0. ( 7 ) 

Eor the foliation-type theories of section IV, it is more conve¬ 
nient to work with covariant Uc, in which case the expression 
for the Einstein equation is 

+ ( 8 ) 

In some cases the Einstein-aether action is supplemented with 
a Lagrange multiplier term enforcing the unit constraint. 


Sx 


d^X X {^gabU°‘U^ — l) . 


( 9 ) 


This will contribute a A-dependent term to the Einstein equa¬ 
tion as well as the aether variations. We will write such terms 
explicitly when they appear; the quantities £ab, d^ab and £a 
are always defined by (4) and (5), with the action S given by 
( 1 ). 

Most of the theories considered in this paper involve scalar 
fields that determine the background structures on which the 
theory is based. We will often make use of the fact that the 
scalar equations of motion are implied by the other field equa¬ 
tions. The proof of this is straightforward: consider an action 
S[gab,Xj 4**] that is a diffeomorphism invariant functional of 
the metric, other tensor fields %, and scalar fields $*. Under a 
diffeomorphism generated by the action varies by 

= (10) 

This variation must vanish for all vectors and since the 
first two terms are zero when the metric and x field equations 
hold, we find that 

= ( 11 ) 


As long as the gradients are (non-vanishing and) linearly 
independent, which can hold for up to four scalars, this im¬ 
plies that the scalar field equations hold, 6S/S^^ = 0. If their 
equations of motion are automatic in this way, we can fix the 
scalars at the level of the action without losing any dynamical 
information. The gauge-fixed action is no longer invariant un¬ 
der the full diffeomorphism group. In our application, where 
a vector field u“ is constructed using scalars, the gradients 
of the scalars must be linearly independent in order for to 
be nonsingular and non-vanishing. Thus, for the physically 
relevant configurations, those scalar field equations are auto¬ 
matic.^ A scalar Lagrange multiplier, on the other hand, need 
not have non-vanishing gradient, so its field equation (i.e. the 
corresponding constraint) should be imposed directly. 

Einally, we recall that in any diffeomorphism invariant the¬ 
ory, some of the field equations are constraints on initial data, 
rather than evolution equations. Eor Einstein-aether theory, 
which contains a dynamical, contravariant vector field the 
quantities 

( 12 ) 

contain no more than first partial derivatives with respect to 
t, for any choice of t and the remaining three coordinates 
[14, 15]. When the field equations are satisfied, we have 
= —Xu°‘Ub and £b = 2Xub, where A terms arise from 
the Lagrange multiplier term (9). The A terms thus cancel, 
so in Einstein-aether theory vanishes on shell. When t 
is a time evolution coordinate, these constraints thus restrict 
the allowed initial data. Eor a covariant aether vector Ua, the 
expression for the constraint has a different appearance, 

=yat{2F\-e^Ub), (13) 

but is in fact the same as a consequence of (6). 

More generally, in the various theories we consider here, 
although the quantities (12) (or (13)) will have no higher than 
first f-derivatives (in appropriately adapted gauges), some or 
all of them may not vanish when the field equations hold, be¬ 
cause the u°‘ (or Ua) field equation per se is not imposed. Eor 
this reason, we refer to them generally as “constraint quan¬ 
tities,” rather than as “constraints.” This failure of constraint 
equations to hold corresponds to the lack of diffeomorphism 
invariance of the gauge-fixed action. We shall analyze the 
form of the constraint violation in each case as a means of 
characterizing the extra freedom available in solutions to these 
theories. 


III. TWISTING AETHER: THREADING THEORIES 

A twisting aether flow does not determine a preferred foli¬ 
ation of spacetime by spacelike hypersurfaces, but it does de- 


^ It can happen that the gradients fail to be independent on a set of mea¬ 
sure zero, e.g. on a codimension one surface, but with remaining well- 
defined in the limit as that surface is approached. In that case presumably 
continuity implies that the scalar field equations also hold directly on that 
surface, at least provided the fields are all nonsingulai'. 
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fine a preferred threading of spacetime. In Einstein-aether the¬ 
ory, this threading is specified by an independently varied vec¬ 
tor field modulo local scale, from which the aether 4-velocity 
(3) is constructed with the use of the metric (alternatively, 
one can use a Lagrange multiplier term to enforce the unit con¬ 
straint on u“). In this section we consider three other ways of 
constructing In the first subsection, the threading is de¬ 
termined by three scalar fields which are Lagrangian (comov¬ 
ing) coordinates for the aether. In the following two subsec¬ 
tions, an additional scalar field ip is introduced into both the 
line field and the Lagrangian coordinate constructions of 
The field ip is an independent “clock” that marks time along 
the threads, and is constrained to agree with proper time by 
a Lagrange multiplier term. These constructions are all very 
closely related to each other, but they yield theories that dif¬ 
fer insofar as different integration constants are required to 
determine a solution, corresponding to different initial value 
constraints that are violated. 


A. Fixed threading theory 

A threading can be specified as the curves along which 
three scalar fields (p\ I = 1,2, 3 are all constant. If the the¬ 
ory is to depend only on these curves as one dimensional sub¬ 
manifolds, and not on any parameterization, the action must 
be invariant under all smooth invertible field redefinitions of 
the scalars. 


(14) 

This can be achieved by restricting the action to depend on 
only via the unit aether 4-velocity 

u%p^,gab) = r/\Jl (15) 

with J“ the metric-independent vector density 

ja ~abcd, A ,^3 

■J = e (16) 

where is the alternating symbol, i.e. the Levi-Civita ten¬ 
sor density of weight 1. The vector field defined in (15) is 
invariant under the ‘V-diffeos” (14), since both the numera¬ 
tor and denominator are rescaled by the Jacobian determinant 
det {d(p^/dAY The corresponding action (1) is then a func¬ 
tional of the metric and the three scalar fields. 

Note that the action is quadratic in second derivatives of 
the scalars A ■ This implies that the field equations will be 
fourth order in derivatives of the scalars, and third order in 
derivatives of the metric (arising from the Christoffel connec¬ 
tion terms). However, as explained in section II, we may treat 
the scalars as fixed, not varied in the action, without chang¬ 
ing the dynamical content of the theory. Since the A define 
a threading, we call this the threading theory. In the 
co-moving gauge, where A equal to the spatial coordi¬ 
nates, contains no derivatives, and the field equations aris¬ 
ing from metric variations are of second order. 


1. Relation to zero temperature perfect fluid 

The dynamics of perfect fluids was formulated long ago in 
terms of three Euler potentials A [16, 17], a formulation that 
has recently been fruitfully exploited with the application of 
ideas from effective field theory (see e.g. [18]). In that setting, 
the vector density J“ represents the conserved entropy current 
and I J| is the entropy density in the fluid rest frame. The en¬ 
tropy current is invariant under t/j-diffeomorphisms with unit 
Jacobian determinant. Unlike for our aethereal application, 
full (^-diffeo invariance is not imposed, because the entropy 
density is physically meaningful. The presence of conserved 
particle number necessitates an additional scalar field with a 
shift symmetry in the action. Our “clock field” ip introduced 
below [see e.g. (46)] is directly analogous to this, although 
the corresponding chemical potential u°‘Vatp is required by 
the unit norm constraint to be everywhere equal to unity. 

Lor fluids without conserved particle number, the action at 
first order in derivatives is the integral of minus the energy 
density expressed as a function p{b) of the entropy density 
scalar b = \J\/y^—g. The function p(b) determines the equa¬ 
tion of state of the fluid. The aether fluid has the property that 
at first derivative order its energy density is independent of 
the entropy density, as required by the full (/3-diffeo symme¬ 
try. This happens for a thermal fluid only at zero temperature, 
hence the aether can be considered a zero temperature fluid. 
The action for such a fluid is just proportional to the spacetime 
volume, so the stress-energy tensor at this derivative order is 
nothing but a (possibly vanishing) cosmological constant, mo¬ 
tivating the name “vacuum fluid” for the aether. The dynamics 
of the vacuum fluid is governed at lowest derivative order by 
the action (1) involving the “strain” of the fluid. 


2. Comparison with Einstein-aether theory 

Under variations of the fields, in u°-{A^,gmn) (3) and 
u°'{Aj 9mn) (15), the variation of has both a parallel and 
a perpendicular part. 


Su'^ = SA6u°p. (17) 

The metric variation generates only Su^ while the A™ and A 
variations generate only Su°[. The metric-induced variation in 
both cases is just what is needed to keep a unit vector: 

0 = 5{gmnU^uA = {5gmn)u^u^ + ‘^9maU^5u°' implies 
that hence 

(18) 

The perpendicular part of the variation is given in Einstein- 
aether theory by 

5ul = (C - u‘^Um)^A^ (19) 


and in the fixed threading theory by 
1 


5u\ = 


2|J| 


(c - u^^uA-^^^AijKAbAc (20) 
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Thus the metric equation of motion (7) is the same in terms 
of gmn and u™ in the background threading theory as it is in 
Einstein-aether theory, 

Sab ScUaUl) — 0 . ( 21 ) 

The remaining equations of motion arise in both theories from 
the variation Su^^, and here a discrepancy arises. 

In Einstein-aether theory the equation of motion arising 
from perpendicular aether variation (19) is 

Sj^ = (<5“ - u^u^)Sa = 0, (22) 

while in background threading theory it is 

( 23 ) 

Diffeomorphism invariance of the action implies that the 
scalar equations (23) hold as a consequence of the Einstein 
equation [cf. discussion around Eq. (11)], so they add no 
new information. On the other hand, the perpendicular aether 
equation (22) adds restrictions in Einstein-aether theory. 

To discover the precise relation between the equations of 
motion (22) and (23), note that since u°‘{(p\gmn) is con¬ 
structed covariantly, diffeo variations of its arguments induce 
its diffeo variation = £^u°‘ as a vector field. In particular, 
the perpendicular component of £^u°‘ is equal to the variation 
induced via £^ip^, so the corresponding contributions to the 
variation of the action S[gab, w“((/?^, (?mn)] are also equal, 

J £i£iu'^ = J £K£ig^^. (24) 

Now using £^u°‘ = — and integrating by parts we ob¬ 
tain 

Ji£a£i-£K^%)r=0 (25) 

for all vector fields which yields the identity 

£u£^=£k^%. (26) 

Thus (22) implies (23) (provided again that the gradients (p^ 
are linearly independent), but (23) implies only that the Lie 
derivative of (22) holds. 

We thus see that in the fixed threading theory the Einstein 
equation implies 

£u£i = 0. (27) 

Put differently, instead of the aether equation (22) one has 

= -Ai, (28) 

where the “source term” is a covector density that satisfies 
= 0 and is conserved along the aether flow, 

£upi = 0. (29) 


A transparent way to express the conservation law (29) is to 
use adapted coordinates, , and to choose = t with 

Q = 1, so that the components of the aether 4-velocity 
are all constant, = i5“. Then the components of the Lie 
derivative are just the partial derivatives with respect to r, and 
(29) takes the simple form 

drPi = 0 . ( 30 ) 

The three components pf are then just constants of integra¬ 
tion on each thread, while p^ vanishes identically. The free¬ 
dom to choose these integration constants different from zero 
is what distinguishes the fixed threading theory from Einstein- 
aether theory. Eor lack of a better name, we shall call p:^ the 
vector source density (VSD). 

The identity (26) shows that the equation of motion is 
“weaker” than the perpendicular it™ equation of motion, but 
this discrepancy remains a bit mysterious, since it would seem 
that variations of produce all possible perpendicular vari¬ 
ations of u"*. Of course the difference must arise because 
occurs in the action with an extra derivative, but why exactly 
is that important? The answer lies in the boundary conditions. 
When we drop boundary terms we are holding fixed at the 
boundaries, in particular the initial and final boundary. This 
entails an integral constraint on the variations (the end¬ 
points of each thread are fixed), which translates into the fact 
that the p^ variations imply only the time derivative of the 
equation of motion. 

A simple example serves to illustrate this point. Suppose 
we have a mechanical system in one dimension with La- 
grangian L{x,x), and we make the replacement x = y, and 
treat y{t) as the basic dynamical variable. Then the action 
variation is 5S = J{SS/Sx)6y = — J{d/dt){SS/Sx)Sy, so 
the y equation of motion is the time derivative of the x equa¬ 
tion of motion. It is weaker than the x equation of motion 
because not all x variations are included in the y-version of 
Hamilton’s principle. Since the initial and final values yip are 
fixed in the y-variations, there is an implicit constraint on the 
integral of x due to the fact that J xdt = J y dt = y 2 — yi- 
We could include this constraint directly in the a;-version of 
Hamilton’s principle with the addition of a Lagrange multi¬ 
plier term A(/ xdt — Ay). The result would be the equation 
of motion SS/Sx = A, where the Lagrange multiplier A is an 
undetermined constant corresponding to a constant external 
force. In the y equation, A corresponds to the extra integration 
constant needed to specify a solution. 

How does a nonzero VSD for the aether field equation 
change the aether theory? In another paper, we find that it 
does not alter the Newtonian limit or static, spherical stars 
(assuming no radial aether component) [19], and by symme¬ 
try homogeneous, isotropic cosmology is unaltered. However, 
it acts as an external force for wave modes, shifting the equi¬ 
librium amplitude away from zero. In the next subsection we 
show that, more generally, the source density integration con¬ 
stants characterize a violation of the initial value constraints of 
Einstein-aether theory. In the following subsection we show 
that magnitude of the source density is diluted as the aether 
expands with the universe, which suppresses its observable 
consequences. 
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3. Initial value constraint violation 

The VSD in (28) suggests that the fixed threading the¬ 
ory requires more initial data than Einstein aether theory, since 
is a freely specifiable initial source for the aether equation. 
This new freedom can be characterized in terms of violated 
Einstein-aether initial value constraints. 

Eor an arbitrary fourth coordinate x° , the constraint quanti¬ 
ties (12) take the form 

= (31) 

When the metric field equation (21) holds, (28) implies that 
these quantities are nonvanishing and instead satisfy 

= -uVa- (32) 

The u“ component constraint u'^Ca^ = 0 holds, since 
= 0’ but the three “perpendicular constraints” are vio¬ 
lated. 

In adapted coordinates, = ip^, the constraint violation 
is preserved in time. This is easiest to see with the choice 
= r, with T the proper time along the threads. Then we 
have = 1, and (30) shows that the components of the 

constraint C^'‘ are preserved in t. In fact the same result 
holds for any choice of the fourth coordinate the condi¬ 
tion = 0 implies that under a change from r to a;°, the 

components of the covector density p^ change only by the 
Jacobian factor dr/dx^, while vP = {dx^/dT)u^. Therefore 
the components of the a;°-constraint (32) in {x^,x^) coordi¬ 
nates are the same as those of the r-constraint in (r, x^) co¬ 
ordinates. Since u'^do = dr, the previous result implies that 
also"^ 

= 0. (33) 

This equation shows that the new freedom takes the form of 
an infinite collection of conserved quantities. The constraint 
violation may be freely specified at an initial time, but remains 
constant at all subsequent times. 

The vanishing of the constraint quantities in Einstein-aether 
theory is a consequence of full spacetime diffeomorphism 
symmetry. The fixed threading theory respects only the thread 
preserving diffeomorphisms, which in adapted coordinates 
take the form t i—>■ t{t,x^), x^ i—>■ x^{x'’). Intuitively, since 
we cannot perform arbitrary gauge transformations of the spa¬ 
tial coordinates as we evolve in time, there should be no 
constraints associated with those diffeomorphisms imposed 
on the dynamics. This is why we find that the spatial 
constraint quantities for evolution along the threads, are non¬ 
vanishing. 

By contrast, for evolution with respect to a parameter that is 
constant on the threads, say x^, all constraint quantities van¬ 
ish, since u°'VaX^ = 0, so the number of initial value con¬ 
straints remains equal to four. This might be expected since 


A covariant version of this argument uses (29) and the identity 
“^u/(u.dxO)[(“ ■ which holds in view of the unit 

density weight and u^flp = 0. 


as we evolve in x^, we can perform both time and spatial dif¬ 
feomorphisms. (That these are required to preserve the fibers 
evidently does not cause the constraints to be lost.) This gauge 
symmetry means the dynamics cannot be fully deterministic, 
so that some field equations must be constraints. 


4. Cosmological evolution of source density 

In homogeneous isotropic symmetry, the VSD necessarily 
vanishes, and the background threading theory is identical to 
Einstein-aether theory. It is natural to imagine some kind 
of fluctuations around the symmetric configuration however. 
Since the VSD arises as integration constants, its power spec¬ 
trum cannot be derived from the properties of quantum vac¬ 
uum fluctuations. At this point we have identified no princi¬ 
ple to select a primordial spectrum of VSD. What we can say 
however is that the amplitude will decrease as the universe 
expands. 

To characterize the amplitude of the VSD we use the scalar 
quantity 

(34) 

An approximate redshift law for k can be easily obtained 
by using for the homogeneous isotropic metric ds^ = 
dt^ — a{t)^dx^dx^, and neglecting the anisotropic corrections 
to the conservation law (29). Then the coordinates {x\t) 
are adapted to u“, and the conservation law takes the form 
dtp:^ = 0, so (34) yields k cx a“^. The physical effects of the 
VSD therefore decrease like those of radiation as the universe 
expands. 


B. Fixed clock theory: a ponderable aether 

In section III A, we introduced three scalar fields that de¬ 
fined a threading. In the co-moving gauge, these scalars have 
the effect of breaking spatial diffeomorphism symmetry when 
fixed at the level of the action. In this section we consider a 
different theory, in which temporal rather than spatial diffeo¬ 
morphism symmetry is broken. This involves introducing a 
“clock” field Ip that defines a preferred notion of time along 
the aether flow. 

Since the clock field ^ is a scalar, we may again fix r/; to a 
background value at the level of the action. In analogy to the 
fixed threading we expect this fixed clock to lead to a viola¬ 
tion of an initial value constraint and therefore to produce, in 
effect, an additional degree of freedom. The constraint viola¬ 
tion in this case is quite analogous to the “dark matter as an 
integration constant” [13] in projectable Hofava gravity. The 
latter is due to the absence of the local Hamiltonian constraint 
in that theory. That constraint normally arises from the vari¬ 
ation of the lapse function N = but in projectable 

Hofava gravity N = N{t) depends only on t. It is therefore 
not varied independently at each point on a constant t surface, 
so the associated local constraint is not imposed. The covari¬ 
ant construction of an aether with a fixed clock given here 
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yields a similar effect. Unlike in the prqjectable Hofava case, 
however, the “dark matter mass current” is conserved in the 
fixed clock aether theory. 

We start as in Einstein-aether theory with a dynamical vec¬ 
tor field A°' but, rather than defining the aether 4-velocity di¬ 
viding by I A|, we define it by 




Aa 


(35) 


where -ip is the clock field. By construction we have u°">p,a = 
1, so r/) is a parameter on the aether flow compatible with u°‘. 
Note that (35) is unchanged under a thread-dependent shift 
^ Ip + V, with V constant along each thread, A°‘Va'v = 0 
(note this symmetry was called a “chemical shift” in the works 
on effective field theory for fluids [18]). The requirement of 
this symmetry precludes standard kinetic or potential terms 
for Ip and, since u°‘ip^a = 1, a term like {u°‘'ip^af' only adds a 
constant to the action. 

Unlike its Einstein-aether cousin u°'{A'^,gmn) (3), 
u°'{A'm^ Ip) is not a unit vector by construction, so we impose 
the unit constraint by adding a Lagrange multiplier term (9) to 
the action, enforcing the relation 

= QmnA'^A^. (36) 


It seems at first that this could be satisfied either by solving 
a first order ODE for -0 on each thread, or by restricting gmn 
(the condition is independent of the scale of A™ so it can not 
be satisfied by restricting that scale). However, solving (36) 
for ip by integrating along each thread would be inconsistent 
with fixing ip at both endpoints in Hamilton’s principle unless 
further constraints on variations of A°‘ and gab are imposed. 
Instead, it is simplest to view the unit constraint as fixing a 
component of the metric in terms of A™ and ip^m- Since pj 
is a scalar field, its equation of motion is satisfied by virtue 
of the other equations of motion (provided ip^rn. ^ 0). It can 
therefore be considered fixed. We call this fixed clock the¬ 
ory since, when the unit constraint is satisfied, ip marks proper 
time on each thread. 

The variation of (35) is given by 

Su- = - u<^fi,m)SA^. (37) 

The Ip equation of motion thus takes the form of a current 
conservation law. 


which in light of (38) becomes 

Sab + -u‘'£cUaUb + - fiUaUb = 0. (42) 

The aether stress tensor thus picks up the extra dust-like con¬ 
tribution, fiUaUb, which is not present in Einstein-aether the¬ 
ory. Together with the conservation equation (38) this sug¬ 
gests the interpretation of jl as the internal energy density 
of the aether, and motivates the descriptive term ponderable 
aether. Notice that even in the absence of the aether terms in 
the action (1), the Lagrange multiplier term (9) alone suffices 
to introduce the aethereal dust stress tensor pUaUb- 

The sign of p is not fixed, so the aethereal dust can con¬ 
tribute negative energy density to the source in the metric 
field equation. Normally the presence of negative energy 
could produce an instability, since the system could evolve to 
a highly excited state of compensating positive and negative 
energy components. However, as a consequence of the con¬ 
servation equation (38), the spatial integral of p, is conserved 
along any bundle of aether worldlines, precluding instabilities 
of that sort. 

The initial value formulation of the fixed clock theory dif¬ 
fers from that of Einstein-aether theory by a violated con¬ 
straint equation. As explained previously for the fixed thread¬ 
ing theory, the quantity Ca'^ defined in (12) has only first t- 
derivatives of the metric and the aether 4-velocity. Eor the 
background clock theory it has second derivatives, since the 
aether 4-velocity (35) involves the derivative of ip. In the 
adapted coordinate “clock gauge” = ip, however, is 
algebraic. 


u°‘{A>^,iP) = A°^/A°, (43) 

so Ca'^ has only first f-derivatives for any choice of t. When 
the and metric equations (40) and (41) hold, we have 

Sb^ = -('w“f.a)AV'.b- (44) 

The right hand side vanishes when contracted with any vector 
tangent to the constant ip surface, so the presence of nonzero 
jl in the aether equation (40) leads to a single constraint vio¬ 
lation. The additional freedom in the theory is parameterized 
by p. If we choose t = ip as the evolution parameter and use 
the clock gauge, (44) takes the form 


{pvJ^),m = 0, p= py/^ = U°- (2XUa - £a^ , (38) 
The A°‘ equation of motion is 

-^^(C-u>.m)(fa-2Au„) = 0, (39) 

which is equivalent to 

£a - 2XUa = -pip,a- (40) 

In this theory, has no metric dependence, so the Einstein 
equation is 

£ab + XUaUb = 0, (41) 


C^) = -jisl. (45) 

If we further choose spatially adapted coordinates, so that 
= bg , the right hand side of equation (45) is constant in a:° 
as a consequence of (38). As in the fixed threading theory, the 
new degrees of freedom are, in this sense, “totally integrable.” 

We have found that for evolution with respect to any coor¬ 
dinate t such that u ■ dt 0, the 0-component of the con¬ 
straint quantities does not vanish. This is because the fixed 
clock breaks time diffeomorphism symmetry. The clock shift 
symmetry remains, but it allows for only a single, time in¬ 
dependent shift, so the ip surfaces cannot be deformed as we 
evolve along the threads. The components of the constraints 
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in the directions tangent to the ip surfaces are preserved, since 
the threads are determined by a dynamical vector field, rather 
than by a background structure. 

The contribution ^UaUb has the form of a pressureless fluid 
source in the Einstein equation, but its divergence is not zero 
when the aether is not geodesic. In homogeneous isotropic 
cosmology, however, it does behave exactly as pressureless 
dust, with /r cx a~^. During an inflationary period fx would 
be exponentially suppressed, so in the standard inflationary 
cosmological model it would presumably be too small today 
to have any observable effect. If there were some way to 
transcend the classical conservation law for /i and generate 
a nonzero value around the time of matter radiation equality, 
it could play the role of the homogeneous dark matter in a 
ACDM model. ^ This leads to the question of how it would 
behave as structure forms. Its nongeodesic character suggests 
that it would not form structure in the manner of geodesic dark 
matter. In another paper [19] we have examined the growth of 
linearized perturbations, and found that if /r were to comprise 
the homogeneous dark matter density at early times it would 
lead to an unacceptably high growth rate on super-horizon 
scales and no growth on sub-horizon scales. 

C. Fixed aether theory 

In the previous two sections, we considered a theory with 
broken spatial diffeomorphisms, the fixed threading, and one 
with broken temporal diffeomorphisms, the fixed clock. In 
this section we combine these features and consider an aether 
theory with broken spatial and temporal diffeomorphisms. 

We now define the aether 4-velocity by 

where J“ is the vector density constructed from defined in 
equation (16), and ip is a scalar clock field. Like the aether 4- 
velocity of the background threading theory (15), u°'{(p^,ip) 
is unchanged under all (/3-diffeos ( 14) and, as in the fixed clock 
theory, it has the clock shift symmetry, ip i—>■ ip+v{(p^) (again, 
this corresponds to the “chemical shift” in the effective theory 
of fluid dynamics [18]). Also, as in the latter theory, it is not 
normalized by construction, so the Lagrange multiplier term 
(9) is again used to impose the unit constraint. 

Since u°‘{(p\ip) is constructed entirely from four scalar 
fields, its variations arise solely from variations of the four 
scalars. As explained above, provided the four gradients 
and ip^a are linearly independent, which is required for them 
to define a threading with parameter ip, the equations of mo¬ 
tion for the four scalars will follow from the Einstein equa¬ 
tion. They can therefore be held fixed in the action. This fixed 


^ If it were generated earlier, e.g. at reheating, it would quickly dominate un¬ 
less fine-tuned to an extremely small value relative to the radiation energy 
density. Mechanisms for generating dark matter after inflation have been 
proposed for the related projectable Hofava gravity [13] and mimetic dark 
matter [20, 21 ] theories. 


aether theory is thus equivalent to Einstein-aether theory with 
an aether vector that is not varied in the action; that is, 
the aether field equation Ea of Einstein-aether theory is not 
imposed. 

When the gab and A equations of motion hold, however, 
the (vanishing) variation of the action S'[pa6jM“,A] with re¬ 
spect to a diffeomorphism generated by f “ is given by J {£a — 
2Xua)£^u°' = J — 2Xua)- Since this vanishes for 

all we infer that 

£ui£a - 2Xua) = 0. (47) 

This is similar to (27) in the fixed threading theory, but in¬ 
cludes the component of £a along Ua- Thus although the 
aether held equation is not imposed, it holds with the addi¬ 
tion of an undetermined, “constant” source term, 

£a 2XUa — fiai £ufia — 0- (48) 

Equation (48) implies 

A = i (m + u^-Ea ), (49) 

with 

Pt = h-'/^=U°‘fia- ( 50 ) 

Therefore, as in the fixed clock theory, the aether stress tensor 
contribution {X/y^—g)uaUb picks up the extra term gUaUb not 
present in Einstein-aether theory. Also, (48) and £uU^ = 0 
imply £ufi = (/iu“),Q = 0, so that p acts like a “dark matter” 
source of gravity that can be interpreted as the internal energy 
density of a ponderable aether. 

When we work in co-moving, clock gauge (x^ = ,x^ = 

Ip), the diffeomorphism symmetry is broken down to time in¬ 
dependent transformations x^ —?> /^{x'’) and —>■ -f 

f{x'^), so we should expect all four constraints to be vio¬ 
lated. When the metric and A equations are satisfied, the x°- 
constraint quantity (12) for the fixed aether theory in these 
coordinates takes the form 

Cio) = 2fO„+uX = -/ic. (51) 

This indeed confirms that all four initial value constraints of 
Einstein-aether theory are violated. This is as expected, since 
for evolution with respect to any parameter that advances 
along the threads, there remains no diffeomorphism freedom 
that would make the dynamics underdetermined. Eor evolu¬ 
tion with respect to a parameter s that is constant along the 
threads, we again And that all constraints vanish. In addition to 
the s-dependent thread preserving diffeomorphisms, the clock 
field’s shift symmetry allows for s-dependent changes in ip. 
Thus, we And four additional initial value freedoms per spa¬ 
tial point, which again by (48) are “totally integrable.” 

IV. FOLIATION THEORIES 

We now turn to theories involving a foliation of spacetime 
by spacelike hypersurfaces. These theories are distinct from 
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threading theories because the aether vector, constructed as 
the unit normal to the foliation, is necessarily twist-free. The 
simplest foliation type theory is khronometric theory, the low 
energy limit of nonprojectable Hofava gravity. After review¬ 
ing its construction in section IV A, we proceed in section 
IV B to add a fixed clock to the theory as was done for the 
threading theories in section IIIB. As in that case, the result¬ 
ing theory exhibits constraint violation, and similarly contains 
a “dark matter” component in the Einstein equation. Finally, 
in section IV C we consider the case where the foliation and 
the clock field coincide. This results in the projectable ver¬ 
sion of Hofava gravity, and we discuss the relation between 
the constraint violation and the “dark matter as an integration 
constant” of that theory [13]. 

A. Khronometric theory 

A twist-free aether can be described by a scalar field T, 
dubbed the “khronon,” whose level sets define the hypersur¬ 
faces orthogonal to the aether 4-velocity [3-5]. In order that 
the theory depend only on the foliation by hypersurfaces, and 
not the values of T, the the action must be invariant under 
monotonic reparametrizations 

T f(T). (52) 

The gradient of T transforms as Ta = {dT/dT)T^a, so the 
numerator and denominator of the aether 4-velocity (2) both 
acquire a factor dT/dT, and these factors cancel. There¬ 
fore the action S[gab,u°‘{T, gmn)] (1) is invariant under T 
reparametrizations. 

Just as in the threading theory, this action is quadratic in 
second derivatives of T, so when T is varied it yields equa¬ 
tions of motion that are fourth order in derivatives of T and 
third order in derivatives of the metric. Again, as explained 
in section II, we may fix T at the level of the action without 
changing the dynamics. In the adapted gauge where T is iden¬ 
tified with one of the spacetime coordinates, we have from (2) 

u“(r,p,fc)=5“^/V^. (53) 

Since the aether 4-velocity is an algebraic function of the 
metric components in this gauge, and the action (1) pro¬ 
duces a second order field equation for the metric (terms with 
more than two spatial derivatives occur in the full Hofava- 
Lifshitz theory). In this formulation, which is equivalent to 
Hofava’s original one [6], the action is invariant only un¬ 
der T -foliation preserving diffeomorphisms, together with T- 
reparametrizations (52). 

We can now examine the constraints for this theory. 
The metric dependence of Uc induces a variation 6uc = 
— ^UcUaUb5g°‘^, so the metric field equation (8) reads 

^ah 'O'c^a^b — 0 - ( 54 ) 

When this equation is satisfied, the constraint (13) is equal to 

pW _ 

'-6 — 


If we choose t = T then, since dT oc u, the right hand 

(T) 

side vanishes, so all the constraints ’ vanish in the adapted 
gauge. 

This is complementary to the situation described in section 
III A 3 for the threading theory. There we found that for evo¬ 
lution with respect to a parameter constant on the threads, the 
constraint quantities vanish. He we find a similar result; for 
evolution with respect to a parameter s that is constant on 
the foliation, the constraint quantities vanish. We can still 
make s-dependent time reparameterizations and spatial dif¬ 
feomorphisms under this evolution, so we expect constraints 
associated with these gauge transformations. If instead we 
were to consider evolution with respect to a different param¬ 
eter s' that is not constant on the foliation, we would find the 
T-component of the constraint violated. This is because the 
foliation cannot be deformed in an s'-dependent fashion, so 
the theory loses that gauge symmetry and the associated con¬ 
straint. 


B. Fixed clock foliation theory 

We can add a fixed clock to the foliation theory by follow¬ 
ing the method introduced above for the twisting aether. We 
introduce the clock field ip and define the aether 4-velocity 
covector as 

Ua{T,gmn,tp) = ™„j.’ , i (56) 

which is constrained by a Lagrange multiplier term to have 
unit norm. The unit constraint requires that the lapse N = 
be equal to {g‘'^T^a'>P,b)~^, which freezes 
one metric degree of freedom (in the adapted gauge, it fixes 
_ gTTy rpjjg equations for both of the scalars T 
and Ip again follow from the Einstein equation, provided Tq 
and ip^a are linearly independent. Although the gradients dT 
and dip are both timelike, they will generically be independent 
whenever the aether is accelerated, since 

£u{u[aip,b]) = a[ai^,b]- (57) 

Thus, the gauge in which both T and ip are set equal to coor¬ 
dinates will generically be nonsingular. 

The Ip field equation gives a conservation law, correspond¬ 
ing to its shift symmetry ip ^ ip + const., 

= (58) 

Defining, as usual, the scalar density pi = (2Am“ — £‘')ua, we 
now have 

A = i (m + UaSn , (59) 

as in the fixed aether case (49). The metric equation of motion 
receives a contribution from the metric variation in Ua, namely 
Suc = —UcU(^aip,b)5g‘'^■ The metric field equation is then 




(55) 


^ab ^Itatl^b i£ ‘2XlJ, ^Ilc^ia'P^^b) — 0- (60) 
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Rearranging to compare with (54), this becomes 

^ab 'bl'c^a^b ~t“ {^a'^b “t“ ‘^^{a'4^^b)^ — (61) 

where tfjj^ = ijj^b — Ub is the projection of ijj^b perpendicular to 
Ub- The p, term gives the difference between this theory and 
the khronometric theory, and it takes the form of non-geodesic 
“dark matter” with momentum density. 

Unlike the previous fixed clock theories, here the clock field 
itself has an effect on the dynamics via the perpendicular com¬ 
ponent of '0 h in the stress tensor (61). Since = 1, V' 

is determined on each thread by its value at one point. Hence 
the value of ip on one spacelike hypersurface must be chosen 
as initial data in order to integrate the equations of motion. 

We now examine the constraints. Using T as the time co¬ 
ordinate, enforcing the Einstein equation (61), and using the 
fact that u cx dT, we find for the constraint quantities 

Cf) = -N-^p-p^b. ( 62 ) 

Thus the T-surface constraint is violated in the ^-component, 
for essentially the same reasons given in section III B for the 
fixed clock theory. The conservation law (58) implies £up = 
0, and we have £ud'ip = 0, so the constraint violation satisfies 

£u = 0. (63) 

This conservation law is more complicated than the analogous 
one for the threading theory (33). In the adapted gauge where 
T and ip are coordinates (and in which the field equations are 
second order in derivatives), the vector cannot be chosen 
to be proportional to dx, since u ■ dip = 1. Instead, we can 
choose a gauge where u = N~^dT + d^, in which case (63) 
becomes 

dTCfi + d^{NCp) = Q. (64) 

Hence, we see that the constraint violation evolves accord¬ 
ing to a first order differential equation that also involves ip- 
derivatives of the metric component N. 

C. Projectable Hofava gravity 

Our final example of a foliation theory is projectable 
Hofava gravity, which can be obtained from khronometric 
gravity (section IV A) by imposing in the action the restriction 
that the lapse function be constant on each foliation surface, 
N = N{T). The aether then takes the form Ua = N(T)T^a- 
An aether of this form satisfying the unit constraint is geodesic 
(u^^VaUb = u^S/bUa = b{u°'Ua) = 0), SO the acceleration 
term CaO^ in the action (1) becomes superfluous. 

The projectability restriction can be implemented by 
adding to the Lagrangian density a Lagrange multiplier term 
^‘^’’N,aT,b- Here is an antisymmetric tensor density, and 
N{g,T) = is the lapse function. In the 

adapted gauge = T, the new constraint term becomes 
with N = ( 5 °°)“^/^. The Hamiltonian constraint. 


which arises from variations of the lapse, now contains an 
additional piece — from the Lagrange multiplier term. 
This term represents a local violation of the Hamiltonian con¬ 
straint. Since the violation is a spatial divergence, the inte¬ 
gral of the Hamiltonian constraint will still be imposed on a 
compact space without boundary. This integrated constraint 
is the generator of global time reparameterizations, and fol¬ 
lows from the global variations of the lapse function N (T) in 
the original 3-1-1 formulation of projectable Hofava gravity. 
With asymptotically flat boundary conditions, the integrated 
Hamiltonian constraint will not vanish, but equal the flux of 
A* *° through the sphere at spatial infinity. This is consistent 
with the fact that the metric is fixed at infinity, so, in particu¬ 
lar, the variational principle does not include global variations 
of the lapse. 

An alternative version of the projectable theory lacks the 
global time reparameterization gauge symmetry, and corre¬ 
spondingly the global Hamiltonian constraint is not imposed 
[22-24]. This version arises from a fixed clock foliation the¬ 
ory (section IV B) when the clock is required in the action to 
be constant on the foliation slices. To implement this we may 
add to the Lagrangian a Lagrange multiplier term X°‘^ip aT^b, 
with an antisymmetric tensor density. The correspond¬ 
ing constraint implies that dip oc dT, so that (56) takes 
the form Ua = ip/\dip\'^■ Having imposed this form for Ua, 
the A“^ constraint can be omitted. At this point T has disap¬ 
peared from the action, its role being taken over by ip£ The 
unit constraint then implies that \dip\ = 1, so that Ua = ip,a, 
and the lapse in the adapted gauge = ip is fixed equal to 1. 
There is no global lapse variation, and no global Hamiltonian 
constraint. This covariant formulation of projectable Hofava 
gravity has been described before in [3, 25].^ 

The Einstein equation in this formulation of projectable 
Hofava gravity results from explicit metric variations alone, 
since Ua is metric-independent, and it reads 

TXUall^b — 0- (65) 

Using the definition of p in (59) we can rewrite this as 

J~ab 2 ^ l^cll'all'b '^pH’al^b — 0 ; ( 66 ) 

which looks like the khronometric theory equation (54) with 
an additional “dark matter” component.* We note a peculiar 
difference from the foliation + clock theory: p appears in (66) 
with the opposite sign as in (61). Thus, in projectable Hofava 
theory, positive p represents negative energy density, whereas 
it gives positive energy density in the foliation + clock theory. 


® We could have kept the role of T explicit, using dtp oc dT to express (56) 
as Ua = T,a/(ip'\dT\'^). The unit constraint then implies il}'\dT\ = 1, 
so that we have Ua = ■>p'{T)T,a- The lapse function is thus given by 
N = N{T) = ip'{T). 

’ For some closely related theories, see [20, 21, 26-28], 

* In the formulation with the integrated constraint, we would find the same 
Einstein equation, with the identification oi p, = 
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The on-shell value of the constraint quantity (13) associated 
with the clock field ip is 


= PiUb, (67) 

so again we find a single constraint violated, due to the pres¬ 
ence of the “dark matter energy density” fi. 

Unlike previous cases considered in this paper, jl is gener- 
ically not conserved along the aether flow in projectable 
Hofava gravity. The conservation equation comes from the 
clock field equation of motion, which reads 




- 2Au" 


= 0 . 


( 68 ) 


Decomposing this equation into parallel and perpendicular 
components, we find evolution equation for jl, 

£ujL = Va£l. (69) 

Since the aether equation of motion is not imposed, this means 
the “dark matter” may be generated or destroyed along the 
flow of Non-conservation of the “dark matter as an inte¬ 
gration constant” was pointed out in [13], where it was sug¬ 
gested that this could provide a mechanism for the generation 
of dark matter during the early universe. In that paper, it was 
assumed that the theory agrees with general relativity in the 
IR, so that the coupling parameters Ci are zero. (Recall that 
the Lagrange multiplier term results in nonzero jl even when 
the aether couplings are zero.) In this limit, £°‘ = 0, so we 
would recover the conservation equation £ujl = 0 were it not 
for the higher derivative terms included in the the full Hofava- 
Lifshitz theory [13]. We note that the non-conservation of jl 
is potentially problematic. Apparently nothing enforces that 
—jl remain positive, so instabilities might arise. ^ 


V. DISCUSSION 

In this paper we studied a variety of aether theories includ¬ 
ing and modifying Einstein-aether theory and the IR limit of 
Hofava-Lifshitz gravity (khronometric gravity), which differ 
only in how the aether is constructed from the independently 
varied fields in the action. When those fields are scalars, their 
equations of motion are implied by the other equations of mo¬ 
tion, so they may be regarded as defining fixed background 
structures. We found that it can be consistent to include such 
background structures, and that they often induce extra de¬ 
grees of freedom owing to the loss of diffeomorphism con¬ 
straints. 

The specific structures we considered were the fixed thread¬ 
ing and fixed foliation, as well as a fixed clock field that could 
be included in either the threading or foliation theories. We 
also considered a non-fixed threading, described by a vector 


® Only when fl is conserved due to the aether parameters being zero has 
it been shown that fluctuations around positive energy backgrounds have 
positive energy [3, 29]. 


field rather than a triple of scalar fields. For the fixed thread¬ 
ing theory, the Einstein equation was unaltered relative to 
Einstein-aether theory, but the perpendicular aether equation 
of motion and the corresponding constraint equations were 
modified by a constant source term. The foliation theory with¬ 
out additional structure is equivalent to nonprojectable Hofava 
gravity. For this theory, the initial value constraints hold in 
the adapted gauge, but the Einstein equation differs from the 
threading-type theories since the aether appears naturally as a 
CO vector Ua ■ 

The addition of the (fixed) clock field ip modifies the Ein¬ 
stein equation by an additional term that has the form of a 
pressureless dust stress tensor, which can be thought of as 
due to an internal energy density of the aether. We therefore 
called such aethers “ponderable.” The fixed clock also leads 
to a violation of the ip component of the initial value con¬ 
straint. When a fixed clock is added to the fixed threading 
theory, we obtain a “fixed aether” theory, which is equivalent 
to describing the aether by a vector field that is not varied in 
the action. When a fixed clock is added to the foliation theory 
and constrained to be constant on the preferred foliation, the 
projectable version of Hofava gravity results. 

The appearance of a dark-matter-like component in the 
Einstein equation is also a feature of the recently proposed 
mimetic dark matter theory [20]. In this theory, the physi¬ 
cal metric gab is constructed from another metric jjab and a 
scalar cp such that the gradient Va</> is unit by construction, 
Qab = c(p^d(p)gab- It was shown in [30] that this theory 

is equivalent to ordinary Einstein gravity, supplemented with a 
scalar field (p that appears in the action only via the constraint 
term imposing that Vc(/) is unit. The discussion of section 
IV C therefore demonstrates that this theory is equivalent to 
projectable Hofava gravity with vanishing aether action, that 
is, with the parameters Ci set to zero. Thus, the mimetic dark 
matter theory can be viewed as a special case of the aethereal 
theories described here. 

The threading theory formalism discussed here resembles 
the Lagrangian description of perfect fluids. The aether, con¬ 
structed from the comoving potentials , acts as a zero tem- 


*** More generally, all the background clock theories discussed here can be 
formulated using the gab metric construction of [20] (see also [29] for 
the case of a vector field), instead of imposing the unit constraint with 
a Lagrange multiplier term. If we define the physical metric as Qah = 
(g'^‘^UaUbj)gab for a covaiiant aether, or as g^b = {gcdu'"u'^)~^gab 
for a contravariant aether, the aether vector automatically has unit norm 
with respect to the physical metric. The action depends only on the 
conformal class of Qab^ so the gab variation gives a trace free equation, 
— {G — T)u^u^ = 0, where and are the vari¬ 
ational derivatives of the action with respect to gab, and the traces G and 
T are their contractions with gab- Foi' the fixed aether theory, there is no 
variation equation, and instead of the usual Einstein equation we have 
only this Einstein equation with an additional source with pressureless dust 
energy-momentum tensor and energy density G — T. Thus we have re¬ 
covered the ponderable aether discussed in the text. The equivalence to a 
mimetic dark matter theory arises in the case when there is no aether ac¬ 
tion. Then is just the matter stress tensor and is conserved when the 
matter satisfies its equation of motion, so the Bianchi identity implies the 
conservation law Va [{G — T)u°"ll^] = 0. The extra term thus behaves in 
this case as geodesic dark matter dust. 






12 


perature “vacuum fluid”, according to the thermodynamic re¬ 
lations developed in [18]. The vanishing of the temperature is 
closely tied with the enhanced symmetry of the aether fluid, 
which includes all (/j^-diffeomorphisms rather than only the 
volume preserving ones. It was mentioned in [31] that such 
an enhanced symmetry is not possible without adding more 
fields, but this is true only for an action that is first order in 
derivatives." The lowest order terms in the aether action in¬ 
volve second derivatives of . A derivative expansion for a 
fluid action was discussed in [18, 32], and the aethereal terms 
invariant under all (^^-diffeomorphisms appear in the latter 
reference. 

The clock held '0 in the aether theories is analogous to the 
phase held introduced in [18] for fluids carrying a conserved 
particle number. In particular, it possesses the same “chemi¬ 
cal shift” symmetry, 0 i—>• 0 -f v{(p^), which, in the aethereal 
case, corresponds to a freedom to shift the initial value of the 
clock along each thread. The scalar y = u°‘Va'4> has a thermo¬ 
dynamic interpretation as the chemical potential for particles 
charged under the shift symmetry. In our aethereal setting, 
however, y is fixed everywhere equal to unity by construction 
of [see, e.g. Eq. (46)], and 0 measures proper time along 
the flow. 

Having elucidated the basic structure of these theories, it is 
interesting to consider the phenomenological consequences of 
the presence of the background structures. In a companion pa¬ 
per [19], we considered some of the astrophysical and cosmo¬ 
logical implications of the source densities for the threading- 
type theories. In particular, we examined whether the new 
component in the Einstein equation for ponderable aethers 
with fixed clocks could play the role of dark matter. The 
two main results of that analysis are that (i) the “aethereal 
dark matter fluid” has pressure, hence does not seed struc¬ 
ture formation on sub-horizon scales, so another dark matter 
component must be present, and (ii) during matter domina¬ 
tion, the presence of a homogeneous ponderable aether en¬ 


ergy density causes problematic growth of the isocurvature 
modes on super-horizon scales. In particular, for isocurvature 
amplitudes of order 10“® at radiation-matter equality (which 
would be the value expected from inflation [33]), the growth 
at large scales becomes inconsistent with CMB and large scale 
structure observations when the ponderable aether contributes 
more than 1% of the homogeneous energy density. On the 
other hand, these results do not apply to the “dark matter as 
an integration constant” [13] in projectable Hofava gravity. 
That theory results from taking the limit Ca —>■ oo,Ci^ oo 
of Einstein-aether theory [4, 9]. In that limit the dark matter 
fluid is pressureless, and the large-scale isocurvature modes 
are decaying. 

Einally, it should perhaps be emphasized that, in a theory 
with conserved “aethereal dark matter” current, the primor¬ 
dial value of the internal energy density of the aether would 
be driven very nearly to zero if there is an early period of 
inflation. This leads to the curious conclusion that the non- 
dynamical and dynamical aether theories could appear to be 
essentially equivalent in their phenomenological predictions. 
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